Square-gradient mechanism of surface scattering 
in quasi- ID rough waveguides 
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We develop a perturbative approach that allows one to study the surface scattering in quasi- 
1D waveguides with rough boundaries. Our approach is based on the construction of an effective 
"bulk" potential of a very complicated structure. The detailed analysis of this potential reveals 
that apart from the well known terms considered in the previous studies, one should keep specific 
terms that depend on the square of the derivative of the boundaries. As was found, in spite of an 
apparent smallness of these square-gradient scattering (SGS) terms, there is a physically important 
region of parameters, in which they can not be neglected. Our approach also demonstrates that 
the contribution of the SGS-mechanism of scattering strongly depends on the type of the rough 
boundaries (uncorrelated, symmetric or antisymmetric). 



PACS numbers: 73.21. Hb, 42.25.Dd 



I. INTRODUCTION 



Wave scattering in surface-corrugated guiding systems 
remains a quite interesting problem, both from the the- 
oretical viewpoint and due to many applications rang- 
ing from the sound and electromagnetic propagation, to 
the conduction of electrons in mesoscopic wires. In spite 
of many important results obtained to date (see, e.g ., 

Refs. p 0, 1 1 i i 0, 1 ©, IE EL M M M EE m) . 

there is a number of open questions to be resolved. One 
of such problems is to understand the interplay between 
different mechanisms of the scattering that are due to 
principally different terms in an effective scattering po- 
tential constructed with the use of various perturbative 
approaches. 

In the recent papers it was shown that apart from 
the well known terms in the potentials appearing in the 
first order perturbation theory, one should keep some 
specific terms that were neglected due to their seem- 
ing smallness. It was demonstrated that these terms, 
although formally belonging to the second order terms, 
can play a decisive role in the scattering in some region 
of parameters. Since these terms are proportional to the 
square of the derivative of surface profiles, they can be as- 
sociated with the square-gradient scattering (SGS) mech- 



In contrast with Ref. 17[ where a particular case of one 
rough surface was considered, below we develop a gen- 
eral approach by concentrating our attention on quasi- 
one dimensional waveguides with two corrugated surfaces 
that are either statistically uncorrelated or correlated. 
Our main interest is in the expression for the attenua- 
tion length L n (also known as the scattering length, or 
the total mean free path), corresponding to specific n-th 
propagating mode. We assume that the total number N 
of these modes (or conducting channels) is large, N ^> 1, 
and that the lower and upper surface profiles (bound- 
aries) of such quasi-lD waveguides are described by two 
random functions £j_ (x) and £j (x) , respectively. 

In what follows, we consider three configurations of a 
particular interest: 

1. The waveguide with the uncorrelated boundaries, 
£j.(x) and £ T (x). 

2. The waveguide with the antisymmetric boundaries, 



(1.1) 
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3. The waveguide with the symmetric boundaries, 

-Z l (x)=b(x)=Z(x). (1.2) 

In our approach we assume that due to the multiple 
scattering of a travelling wave from the rough bound- 
aries, the longitudinal wave number of an n-th propa- 
gating mode can be written as k n + 5k n , where k n is its 
unperturbed value, 



5k n 



In 



i(2L T: 



(1.3) 
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The real part 7„ is responsible for a roughness-induced 
correction to the phase velocity of a given mode. As is 
known, the shift 7 n does not change the transport prop- 
erties of a disordered system. Therefore, our further anal- 
ysis shall be focused on the attenuation length L n only. 

In order to properly take into account the surface ir- 
regularities we use the approach developed in Ref. [T^. 
Specifically, by means of the canonical transformation we 
formally reduce the problem of the surface scattering to 
the bulk one described by the Hamiltonian H , 

H = H {0) +U (1.4) 

with the effective potential U. It should be stressed that 
such a transformation does not mean an equivalence be- 
tween the surface and bulk scattering. The very point is 
that the obtained "bulk" Hamiltonian H is of very spe- 
cific form that does not allow to use standard random 
matrix approaches based on a complete randomness in 
effective potentials. 

The paper is organized as follows. In Sec. [IT] we start 
with the general expressions describing the waveguide 
with any form of rough boundaries. In Sec. IIIII we write 
the Dyson-type equation and derive approximate expres- 
sions for the surface scattering potential in the form, most 
useful for our perturbative approach. In Sec. IIVI we ob- 
tain the fluctuating part of the potential for three con- 
figurations of a waveguide. In Sec. [V] we present a brief 
description of the formalism for obtaining the average 
Green's function developed in Ref. [I7j. Here we also 
specify the conditions of applicability of our method. In 
Sec. lVII we perform a detailed analysis of the attenuation 
length for three specific configurations of the waveguide. 
In Sec. I VIII we report the data that illustrates our ana- 
lytical results and make clear the role of various terms in 
the expression for the scattering potential. In this Section 
we also discuss the so-called "repulsion" phenomenon dis- 
covered in Ref. UH, in view of its relation to the square- 
gradient scattering mechanism. A brief comparison of the 
three configurations, as well as our concluding remarks, 
are presented in Sec. IVIIII 



II. THE WAVEGUIDE WITH GENERAL 
ROUGH BOUNDARIES 

In this section we derive the general expressions de- 
scribing quasi-lD waveguides with any type of the bound- 
aries. Then, we present them in the form, most suitable 
for our analysis of three particular waveguide configura- 
tions. 

Specifically, we consider the standard model of an open 
waveguide of the average width d, stretched along the x 
axis. The lower and upper surfaces of the waveguide 
are assumed to be described, respectively, by the profiles 
z = <j£i(x), and z = d+a^(x). Here a is the root-mean- 
square roughness height that is assumed to be identical 



for both boundaries. In other words, the waveguide oc- 
cupies the region 

— oo < x < oo, cr£i(x) < z < d + a£-\(x) (2.1) 

of the (x, z)-plane. The fluctuating width w{x) is defined 
by 

w(x)=d+*[Z 1 (x)-Z l (x)}, (w(x))=d, (2.2) 

where the random functions d{x) (with i =1, f ) describe 
the roughness of the corresponding boundary. 

In what follows we consider the model for which both 
of the boundaries are statistically identical. They are 
assumed to be described by a statistically homogeneous 
and isotropic Gaussian random processes with the zero 
mean and unit variance, 

(£<(*)) = 0> (& 2 (z)> = 1, (2.3a) 
(Zi(x)Z i (x'))=W(x-x'). (2.3b) 

Here the angular brackets stand for the statistical av- 
eraging over different realizations of the surface pro- 
files £,i(x). We also assume that the binary correlator 
W(x) decreases on the scale R c , with the normalization, 
W(0) = 1. 

For three configurations we analyze below, the cross- 
correlations are as follows. In the first case of the uncor- 
related boundaries we obviously have, 

(tl(xMx'))=0. (2.4) 

The waveguide with the antisymmetric boundaries im- 
plies: 

(tl(x)^(x'))=W(x-x'). (2.5) 

As for the symmetric boundaries case, the opposite 
boundaries hold the correlation equality: 

(tl{x)^(x'))=-W(x-x'). (2.6) 

The roughness-height power (RHP) spectrum W(k x ) is 
defined by: 

/oo 
dxexp{-ik x x)W(x). (2.7) 
-oo 

Since W(x) is an even function of a;, its Fourier transform 
(|2.7|) is an even, real and non- negative function of k x . 
The RHP spectrum has the maximum at k x = with 
W(0) ~ R c , and decreases on the scale R^ 1 - 

In order to analyze the surface scattering problem we 
shall employ the method of retarded Green's function 
G(x, x'] z, z'). Specifically, we start with the Dirichlet 
boundary-value problem 

(S + ^ + fc2 )^'"' Z ' Z/) 
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= 6(x-x')5{z- z'), (2.8a) 

g(x,x';z = ai L (x),z') = 

G(x,x';z = d + cr^ix),^) = 0. (2.8b) 

Here <5>(a;) and 5(z) are the Dirac delta functions. The 
wave number k is equal to uj/c for an electromagnetic 
wave of the frequency u and TE polarization, propagat- 
ing through a waveguide with perfectly conducting walls. 
For an electron quantum wire, k is the Fermi wave num- 
ber within the isotropic Fermi-liquid model. 



A. Effective surface scattering potential, U 



The equation (|2.8a|) does not contain any scattering 
potential, therefore, the perturbation is hidden in the 
boundary conditions (|2.8b|) . In this case it is convenient 
to perform the canonical transformation to new coordi- 
nates, 



Xnew — X ld-> 

d[z id - cr£,i(xoid)} 

W[x ld) 
d[z id- <j£,l{x ld)] 



d + <j[£,]{Xold) - £l(Xold)Y 



(2.9a) 



(2.9b) 



in which both boundaries are flat. In accordance with 
Eqs. (|2.9[) . the transformations for the first and second 
derivatives with respect to the old coordinates are given 
by (below for convenience we drop the subscript new for 
x and z): 



d 



dx d dz d 



dx . 



■id 



dx, 
d_ 

dx 



old. 



dx dx id dz 

z^ ^rid—z) 



w(x) w(x) 
dx d dz d 



dz Q 



Id 



dz id 
d 



dx 

d 



dz id dz 



w(x) dz 



dx; 



id 



d 


~a&{x) 
w(x) 


<(*) 


dx 


w(x) 


d 


w(x) 


a£l(x) 


dx 


w(x) 



(d- 
(d- 



d 2 



dz; 



hi 



w 2 (x) dz 2 



d_ 

dz ' 

(2.10a) 
(2.10b) 

d_ 

dz 
d_ 

dz 

(2.11a) 
(2.11b) 



Correspondingly, we introduce the canonically conjugate 
Green's function, 



r , i i\ \/w{x)w(x') , / W010 s 

y new (x,x ;z,z) = g o id{x,x ;z,z). (2.12) 

Here the prefactor before g o id is due to the Jacobian of 
the transformation to new variables. As a result, we ar- 
rive at the equivalent boundary-value problem governed 
by the equation, 



U(x, z) g(x, x'; z, z') = S(x - x')6(z - z'), (2.13a) 

(2.13b) 



G(x, x ; z = 0, z ) = 0, 
G(x, x'; z — d, z') = 0. 



Here we also omitted the subscript new for 
Qnew{x,x';z,z'). As a result of the transformation, 
we have formally reduced the problem of surface scat- 
tering to that of the bulk one specified by the effective 
potential U(x, z), 

d 2 1 d 2 

w 2 (x)\ dz 2 



U(x,z) 



1 



w(x) 



t>i(x) + 2t\(x)-^ 



1 d 

2 +Z dz 



z»(x)+2e l ( X )- 



2 + {z - d) Yz 



e; 2 (^) 




3z 



dz 



dz 2 



3 d d 2 

~+3(2z-d)— + 2z(z-d)- 1[ 
Z oz Uz 



(2.14) 

Note that the prime over the function £j(x) denotes a 
derivative with respect to x. It is important to stress 
that Eqs. (|2.13[) and (|2.14[) are exact ones, and valid for 
any form of the surface profiles £i(x) and t;-\(x). 

As is expected, the general expression (12.141) reduces 
to that obtained previously in Ref.[l7j for a particular 
case of one rough boundary, if one substitutes £| (x) = 
and £|(x) = £(x). 



III. DYSON EQUATION 

Here we start with the non-local Dyson equation in 
fez-representation, 

g(k x , k' x ) = 2ir8(k x - k' x ) g (k x ) 
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■9o(k x ) J -^E(k x ,q x )g(q x ,k' x ) 



(3.1) 



with go(k x ) = k z cot(k z d) being the unperturbed pole- 
factor. In this representation the effective surface scat- 
tering potential 5(fc a , k x ) is given by the expression, 



E(k x ,k' x ) = / dxi 



dz\ exp (— ik x x\) — — - z 1 ^ 



xU(xi, z\) exp (ik x x\) 



sm(k z zi) 



(3.2) 



Here k x and k z are the lengthwise and transverse wave 
numbers, 



k z — k z {k x ) — yfe^ fc 2 ", fc 2 — k z {k x ). 
Their unperturbed eigenvalues are k n , 
k n = \Jk 2 — {nn/d) 2 1 n = 1, 2, 3, 



(3.4) 



and Trn/d, respectively. The total number A^ of propa- 
gating waveguide modes (conducting channels that have 
real values of k n ) is determined by the integer part [. . .] 
of the ratio kd/ir, 



N d = [kd/n] 



(3.5) 



After the substitution of the expression (|2.14| for U 
in Eq. (|3.2[) . we realize that 'E(k x ,k' x ) consists of three 
groups of terms of a different nature. The first one has 
the factor [1 — d 2 /w 2 {xi)\, the second and third groups 
have the factors a/w(xi) and a 2 /w 2 (xi), respectively. 
One needs to note that while itself the kernel S(fc x , k' x ) is 
Hermitian, the latter two parts are non-Hermitian ones. 
In order to present each part of S(fc x , k x ) in an Hermitian 
form, we perform the integration by parts for the term 
containing £"(x\) and £"(x\) in the second group. As a 
result, we arrive at the final exact form for the perturba- 
tion potential, 



/OO r-d 
dxi / dzi exp {—ik x x\) 
-oo Jo 



sm(k z Zx) 
kz 



w 2 (x\) 



x«t(n) 



|t + [*(k + K)} 

oz{ w(xi) 

d 



w 2 (xi) 



dzi 
1 



8 



4 azi 



1 / « 9 

2 + (2l - d fe 



a 2 



1 5 
OZ\ 



1 9 , <9 2 

4 +2(zi - d te +( ^- d) V 



2^+ 2(221 



+2^1 (zi - d) 



<9V 



exp (ifc^i) 



sm(/c z Zi) 



(3.6) 



This equation has a peculiar structure, very useful for a 
further analysis. The kernel written in this form consists 
also of three groups of terms, however, now they can be 
associated with different scattering mechanisms. 

Since we are interested in the averaged Green's func- 
tion, we have to calculate the binary correlator of 
S(fc x , k' x ). Therefore, in order to avoid very cumbersome 
calculations, it is reasonable to make a simplification of 
E(k x , k' x ) that does not destroy the Hermitian structure 
of each group of terms. Taking also into account that 
our main interest is in the case of small surface corruga- 
tions a <C d, we can do the following: expand the factor 
(3 3) 1~ d 2 /w 2 (xi) sa 2a[£,-\ (xi) — ^i(xi)]/ d in the first term of 
Eq. (|3.6p and to put w(xi) w d in all the others. In such 
a way, we derive the following approximate expression for 
the surface scattering potential, 



z,(k x ,k x )Ri / dxi / dzi exp {—ik x x\) — 

J-oo JO K z 



T teT(zi)-to)]^ + ^(fc*+*4)] 



1 9 



1 / 9 

2+^fe 



- + 2Z!— + z 2 



S«i ' ~ l Or 2 



i + 2(2*-d)A 



, ., , . sm(k'zi) , 
exp(ik' xXl ) — 77-^ (3.7) 



The expression Q3.7[) contains three types of scatter- 
ing terms of a different physical meaning: the terms de- 
pending on the roughness height o~£,i{x\), derivative of 
the roughness height a^ixt), and square of the deriva- 
tive of the roughness height a 2 t;' 2 (xi). In this connection 
one can speak about the amplitude, gradient and square- 
gradient scattering (SGS) mechanisms. Note that these 
terms are separately associated with each of two surfaces. 
Apart from this, there is a specific group of terms that 
are due to the inter-correlations between two surfaces. 
In what follows we refer to these terms depending on the 
product a ^(x)^i(x), as to the gradient inter-product 
scattering (GIS) mechanism. 

It should be stressed that the GIS-terms, due to their 
proportionality to a 2 , were neglected in previous studies 
of the transport properties of surface-disordered waveg- 
uides. A similar situation occurs for the SGS-terms that 
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have been analyzed only recently in Ref. [It], [H|- As 
we show below, all these terms can give the main contri- 
bution in the scattering, although formally they may be 
treated as negligible ones. 

We start our analysis with the presentation of the ker- 
nel as a sum of its average and fluctuating parts, i.e., 



a (k x , k x ) — (a(k x ,k x )) + z,(k x ,k x ). 



(3.8) 



It can be tested that the average, (S(fc x , k' x )), contributes 
only to the real part j n of the complex renormalization 
5k n of the lengthwise wave number k n (see Eq. (|1.3[> ). and 
therefore, does not change static transport properties of 
the surface disordered waveguide. Thus, we will omit 
it because our interest is in the attenuation length L n , 
not in 7„. Therefore, it is necessary to identify the terms 
that have a non-zero mean- value in the general expression 
(EH)- 



IV. FLUCTUATING SURFACE SCATTERING 
POTENTIAL 

A. Waveguide with uncorrelated boundaries 

One can see that the mean value of the kernel in 
Eq. (|3.7p is defined by the SGS-terms that are propor- 
tional to er 2 ^' (x). In order to construct the fluctuat- 
ing part of E(k x ,k' x ), we introduce the square-gradient 
functions Vj(x) for the lower (i =J.) and upper (i =f) 
boundary, by subtracting the mean values, (£j 2 (x)) — 
(^ 2 (x)) = (e' 2 (x)),fromef(x), 

Vi(x) = g(x)-{Z'\x)), (V i (x))=0. (4.1) 

The functions Vj (x) play an important role in our further 
consideration. In accordance with the Gaussian nature 
of the surface-profile functions £i(x), there are no corre- 
lations between Vi{x) and both £i(x) and £-(x), 

^i(x)Vi(x')) = , &(x)Vi(x')) = , (4.2) 

as well as between the product £',(x)^(x) and Vi(x), 

The binary correlators for Vi(x) are given by 
(V^Viix')) = 2(&x)£(x')) 2 = 2W" 2 (x - x'). (4.3) 

Correspondingly, the Fourier transform of Vi{x) reads, 
{V i {k x )Vi{k , J)=A%8{k x + k' x )T{k x ), (4.4) 

due to the standard relation, 

/>OC 

Vi(k x ) = I dxexp(-ik x x)Vi(x). (4.5) 



Here we have introduced the roughness-square-gradient 
power (RSGP) spectrum T(k x ), 

/>oo 

T(k x )= dxcxp(-ik x x)W" 2 (x). (4.6) 



We also introduce the Fourier transform of the product 
£{( x )£.\( x ) an d its binary correlator , 

/oo 

dxexp(-ik x x)£[(x)^(x), (4.7) 
-OO 

(V^V^K)) = 2n6(k x +k' x )T(k x ), (4.8) 
as well as the Fourier transform of £j (x) , 



£i(k x ) 



dx exp(— ik x x) £,i{x). 



(4.9) 



It should be pointed out that by the integration by 
parts the power spectrum of the roughness gradients 
Ci(x) can be reduced to the RHP spectrum W(k x ). How- 
ever it is not possible to do for the RSGP spectrum T(k x ). 
This very fact reflects a highly non-trivial role of the SGS- 
and GIS-terms, since they can compete with the other 
terms, although they are proportional to er 2 . 

We can write now the fluctuating part of the total scat- 
tering potential in the form, 

E(k x ,k' x ) = E^(k x ,k' x ) + E n (k x ,k' x ) 

+ ~2i(kx, K) + E 2i (k x ,k' x ) + E n (k x ,k' x ). (4.10) 

Here the first term is associated with the terms depend- 
ing in Eq. ()3.7() on a£,i(x) and a^(x), therefore, on the 
amplitude and gradient scattering mechanisms, respec- 
tively, 



sin(fc z zi) 





{k x , k x ) — 




/•oo r-d 


(7 


\ dxx \ 


dj 


-oo JO 






x i 






"1 d ' 


X 





ex\){ik' x xi) 



sin(fc^,zi) 



(4.11a) 



Si|(fc„, k nl ) = 

-o-{S nn ' + cos[7r(n - n')](l - <W)Rt - K>). 



(4.11b) 



The second term in Eq. (|4. 10[) corresponding to the lower 
boundary, is of the same form, 



^n{k x ,k' x ) = 



dxi / dz\ exp(— ik x xi) 



sin(fc z Z!) 
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1 / * 9 

2 + {z ^ d) 0^ 



exp^fc^x) sin ( fc ;^) (4.12a) 



(4.12b) 



Next two terms in Eq. (14. 1 0[) are the SGS-terms, with 
Ci 2 ( x ) being replaced by Vi(x), 

S 2 t(^,^) = 



d 2 . 
xV T (xi) 



1 



c9 



^ 2 



sin^Zi) 



(4.13a) 



ct 2 ^ 



1 



3 (27m) 2 



4 n 2 + n A 



7r 2 (n 2 — n' 2 ) 2 



x cos[Tr(n - n')\ (1 - 5 nn >) ^ F T (fc„ - fe n ,). (4.13b) 



x exp^k^xi] 



sin(fc^i) 



(4.15a) 



2(7rn) 2 



8 n 2 + n' : 



7r 2 (n 2 — n' 2 ) 2 



x cos 2 [7r(n - n')/2] (1 - <W) \ V t (k n - k n ,). (4.15b) 



The kernel S(fc x , fc^,) written as the sum (|4.10j) of spe- 
cially designed terms, has a very convenient form. First, 
all the terms are chosen to have zero average. Second, 
since the functions £i(x), £,i(x)^(x) and Vi(x) do not 
correlate with each other, their Fourier transforms are 
also non-correlated: 



(Uk x )v i (k' x )) = o. 
(i i (k x )v l (k' x )) = o. 

(V i (K)V l (k x )) = 0. 



(4.16a) 
(4.16b) 
(4.16c) 



Therefore, all the terms in Eq. (|4.10[) determining the 
potential, are also non-correlated. Their autocorrelators 
are, 

("li {km kn')^n {k n i , k n " )) = 2ir5(k n k n ")Qi (k n , k n ' ), 

(4.17a) 



Ql(k n ,k„') = a 2 W(k n - k„'), 



(4.17b) 



E 2i (k x ,k' x ) = 

xd JO K z 



d 2 . 



d 



d 2 



- + 2( Zl -d) ^— + ( Zl -dy— . 



x exp(i^xi) 
^2l(k n , k n ') : 



sm(k' z Zi) 



dz\ 



(4.14a) 



2 



1 



1 



3 (27m) 2 
4 n 2 + n' 2 



5nn' 



7r 2 (n 2 — n' 2 ) 2 



(l-S nn ')}V l (k n -k n '). (4.14b) 



Finally, the last term in Eq. (|4.10p is due to the GIS 
terms in Eq. (|3.7p . 



E 2 l(kx,k' x ) — 

^2 roo rd sin(fc z zi) 
I dxi dz 1 exp(-ik x xi) - 

» JO fc z 



d 2 

^{(xi)^(xi) 



d 



d 2 



- + 2{2z x - d) — + 2z 1 (z 1 -d)—^ 



dz\ 



dz 2 



(^2i(kn,k„')E 2 i{k n ',k n »)) = 2ir6(k n - k n „)Q 2 (k n ,k n ,), 



<r A d 2 

Q 2 {km k n ') — - T(k n k n i) 



1 1 



3 (27m) 2 



16 (n 2 + n' 2 ) 2 
7T 4 (n 2 — n' 2 ) 4 



{1-Snn') >, 



(4.18a) 



(4.18b) 



(S 2 j(fc n , fc n /)S 2 |(fc„',fc„")) = 2n8{k n - k n ,f)Q 2l (k ni k n >) 

(4.19a) 



a*d 2 

} 2 \(k n ,k n ') — - T(k n k n i) 



3 2(7m) 2 



64 (n 2 + n' 2 ) 2 , 
7T 4 (n 2 - n' 2 ) 4 



ir(n - n')/2] (1 - 5 m >) 



(4.19b) 



Note that in spite of the absence of correlations be- 
tween the upper and lower surface profiles, there 
are the correlations due to their derivatives since 
{S 2t (k n ,k n ')E 2t (kn',kn")) ^ (see Eq. gHH). 
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As a result, the correlator of the fluctuating scattering 
potential (|4.10p in the /^-representation reads, 

fe n /)ii(fc n /, fe n //)) = 27T(5(fc„ — k n n)Q\k n , fen'), 

(4.20a) 

<3(fe„, fe„/) = 2Q 1 (k n , k n >) + 2Q 2 (k n ,k n ,) + Q 2 i(k n , k n >) . 

(4.20b) 

The factor 2 in the first and second terms arises due to 
independent contributions, Qi{k n ,k n >) and Q2(k n ,k n >), 
from the lower and upper boundaries. 



B. Waveguide with antisymmetric rough 
boundaries 

In this configuration both boundaries of the waveg- 
uide are identical, see Eq. (jl.l[) . therefore, antisymmetric 
with respect to the central line z — d/2. For this specific 
case the width of the waveguide (the distance between 
upper and lower boundaries) remains constant along the 
x— direction. Instead of Eq. f|3 . T|) we have to start now 
with the exact expression ( (|3.6[) ) by substituting the re- 
lation ( (|l.lj) ). As a result, one gets the relatively simple 
expression for the scattering potential, 



^(k x ,k x ) = / dxi I dzi exp {-ik x x\) 

J-oo JO fc z 

(4.21) 

One can realize that all terms proportional to a 2 can be 
written as one SGS-term with a non-zero mean value. 
By making use of the functions defined by Eqs. (|4.1|) . 
(|4.5j) and (|4.9|) with i =|, in what follows we omit the 
subscript i. 

Now the fluctuating part of the scattering potential is, 



E[k x ,k' x ) = Ei +E 2 (k x ,k' x ). 



(4.22) 



The first summand related to the gradient scattering 
terms in Eq. (|4.21[) has the form, 

„ . , f°° f d . . sm(k z z{) 

ax(k x ,k x ) = a / dxi I dz% exp(-ifc :E xi) 

J-oo Jo k z 

x [i(k x + k' x )} ?( X1 ) ^ expCi^n) ^fel (4.23a) 

Si (*;„, Aw) = 2crsin 2 [7r(n - n')/2] £(fe„ - fen')- ( 4 - 2 3b) 

The second summand is due to the SGS-term, however, 
with the function V(x\) instead of ^ ,2 {xi). With the use 



of the Fourier transform (|4. 5|) for the operator V(x), one 
gets, 



/oo r-d 
dx\ I dziexp{-ik x xi) 
-oo JO 



- , sin(fc^i) 



d 

xa 2 V(xi) exp(ik' x xi) 



^2{krn k n ') — 2 ^nri' V(jc n fen')' 



sin(fc z zi) 
k z 

(4.24a) 
(4.24b) 



Note that the expressions (I4.23P and (|4.24[) for the scat- 
tering potential are still exact. In contrast with other 
configurations discussed below, the expansion in a/d- 
smallness is not required because of the constant width 
of the waveguide. 

The correlation properties of the scattering potential 
are described by the following formulas. Due to the con- 
dition l|4.16ap . the terms Ei(k x ,k' x ) and E 2 (k x ,k' x ) are 
non-correlated, 

(Ei(fc XI A4)S 2 (A4,A4')>=0. (4.25) 

On the other hand, the autocorrelators are non-zero, 

(iii {k n j k n i (fen' , k n " )) — 2ir8(k n k n ")Qi(k n , k n i) , 

(4.26a) 

Qi(fe„, fev) = 4a 2 W(k n - k n ,) sin 4 Kn - n')/2]. 

(4.26b) 



(3 2 (fc n ,fe n ')S2(fc n ',fe n »)) = 27T(5(fc„ - k n ")Q 2 (k n ,k n '), 

(4.27a) 

<32(fen,fen') = (o- 4 d 2 /2)T(k n -k n ,)6 nn ,. 

(4.27b) 

Finally the correlator of the fluctuating part of the total 
scattering potential S(fe x , k' x ) has the following form, 

(^(fenj k n /)n(k n i, fc n ")) — 2 r K5(k n — k n ii)Q(k ni k n >) 7 

(4.28a) 

Q(fen,fen') =<3l(fen,fen') + Q2(fen,fen')- (4.28b) 

C. Waveguide with symmetric rough boundaries 

Let us now consider the waveguide with symmetric 
boundaries, see Eq. (|1.2[) . therefore, with the width de- 
fined by 



j(x) =d + 2<r£(x), (w(x)) = d. 



(4.29) 



By substituting the symmetric profile into Eq. (|3.7|) . we 
obtain for the surface scattering potential, 



"fei k x ) 



oo fd 



oo JO 



dx\ I dzi exp {—ik x x\ 



sin^zi) 
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f jggfo) d 2 
\ d dz\ 



^ ,2 (*i) 

x exp (ik' x xi) 



1 + 4 (2*i - d) ^- + (2*i - d) 2 



(4.30) 



As for the previous configuration, see Eq. Q4.2ip . the 
SGS-terms corresponding to the lower boundary and the 
GIS-terms are written in terms of £(x) = £f(x). 

The potential in Eq. (|4.30[) can be written as the sum 
of two terms related to substantially different scattering 
mechanisms, 



E(k x ,k' x )=E 1 {k x ,k' x ) + E 2 {k x ,k' x ). 



(4.31) 



Here the first summand is related to the terms propor- 
tional to a in Eq. (|4.30p (the amplitude and gradient 
scattering terms), 

„ . f°° f , . sin(fc z zi) 

a 1 [k x ,k x )= / dxi / dzi exp(-ik x xi) 

J — oo Jo k z 



1 + (2zi - d) 



d 



exp(ik' x xi) — -tj-^- (4.32a) 



Hi(fc„, fe„/) = -2cr{5 nn / + cos 2 [7r(n - n')/2] 

X(l - <W)}£( fc n ~ 



(4.32b) 



The second summand is associated with the fluctuating 
part of the SGS-terms in Eq. (|4.30p . and consequently, 
with the function £' (x) being replaced by V(x), 

- „ /"^ , /" d , / - s sin(fc zZl ) 
^2(k x ,k x ) = - / ctai / etei exp(-zfe x a;i) 



d 2 



1 + 4 (2 Zl - d) ^- + ( 2zi - d) 2 ^ 



dz 2 



x enp(ik x xi) 



^(km fen') 



n 2 + n' 2 
(n 2 — n' 2 ) 2 



sm(k' z Zi) 



(4.33a) 



<r 2 d 



1G 



cos 2 [^(n - n')/2](l - <W) - 

(4.33b) 



The correlation properties of the introduced surface 
scattering potentials Si(fc„,fc„/) and S2 (fen, fen') are de- 



scribed by the correlator, 

("l [kn ? fen' [ten' 7 k n " )) = 1iT8(k n k n ")Qi(k n , k n >), 

(4.34a) 

<3l(fe„,fc n ') = 4cr 2 VK(fc„ - fc„/) 

x{<W + cos 4 [^(n - n')/2](l - <W)}- (4.34b) 



(S 2 (fc„, fcn')S 2 (fc„', fen")) = 27T(5(fc„ - fc n " )Q 2 (fc n) fen')) 

(4.35a) 



cr 4 d 2 

?2(feri7 fe«') = ~ T(k n — k x ) 



(ttti) 2 



+ 



256 (n 2 + n' 2 ) 2 



7r 4 (n 2 — n' 2 ) 4 



cosV(n-n')/2](l-<W)j. 



(4.35b) 



Thus, we arrive at the expression, 



(£2(fe n , fen')s(fe n ' > fe«")) — 27T(5(fc n k n ii)Q\k n ,k n ')f 

(4.36a) 

Q(k n , fc n /) = Qi(k n , k n i ) + Qi(k n , k n '). (4.36b) 



V. AVERAGE GREEN'S FUNCTION AND 
ATTENUATION LENGTH 



Now we replace the problem for the random Green's 
function g(k x ,k' x ) with the problem for the Green's 
function (g(k x ,k' x )) averaged over the surface disorder. 
To perform the averaging of Eq. (|3.1[) with S(fea;,fe^), 
given for the uncorrelated, antisymmetric and symmetric 
boundary configurations by Eq. (|4.10[) . Q4.22[) and Q4.31[) . 
respectively, and obtain (g(k x ,k x )), we can apply one of 
the standard and well known perturbative methods. For 
example, it can be the diagrammatic approach developed 
for surface disordered systems [l|] , as well as the technique 
developed in Ref. d. Both of the methods take adequately 
into account the multiple scattering from the corrugated 
boundaries, and allow one to develop the consistent per- 
turbative approach with respect to the scattering poten- 
tial. After quite cumbersome calculations (see Ref. [ItJ 
for details), we obtain the average Green's function in 
the normal-mode representation, 

N d / l\ 

(G(x, x'; z, z')) = ^2 sin (~~r) sin ( ) 

n— 1 ^ ^ 



exp(ifc„|a; — x'\) 
X — —ri ex P 

1 Kifi (X 



\x — X 
2L n 



(5.1) 



The quantity L n is the wave attenuation length or elec- 
tron total mean-free-path of the nth conducting mode. 
This quantity is determined by the imaginary part of 
the self-energy and describes the scattering from the nth 
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mode into other propagating modes, 



1 

La 



(irn/d) 2 
k 7l d 



g (nn'/d) 2 

n' = l 



k 7} /d 



x[Q(k ni -k n ,) + Q(k n ,+k n ,)}. (5.2) 

Here the explicit expressions for Q(k n ,—k n >) and 
Q(k n ,+k n >) are given by Eqs. (I4.20b|) . (14.28bj) and 
(|4.36bp for the waveguide with uncorrelated, antisym- 
metric and symmetric boundaries. 

In deriving Eqs. (|5.ip and (|5.2p we have made the 
following simplifications. First, the self-energy in the 
Dyson-type equation for the average Green's function 
was obtained within the second-order approximation in 
the perturbation potential. In terms of the diagrammatic 
technique this is similar to the "simple vortex" or, the 
same, Bourret approximation [20] which contains the bi- 
nary correlator Q(k x ,k' x ) of the surface scattering poten- 
tial and the unperturbed pole factor go(q x ). Second, in 
order to extract the inverse attenuation length from the 
self-energy, we substituted the lengthwise wave number 
k x by its unperturbed value k n . The justification of this 
substitution, as well as the conditions of applicability re- 
sulting from the above simplifications, were discussed in 
Ref. [17[. Specifically, it was shown that the domain of 
applicability is restricted by two independent criteria of 
the weak surface scattering, 

A n = k n d/(im/d) < 2L n , R c < 2L n , (5.3) 

Here A n is the distance between two successive reflections 
of the nth mode from one boundary to the opposite one. 
The conditions in (|5-3[) imply that the wave is weakly at- 
tenuated on both the correlation length R c and the cycle 
length A„. They implicitly include the requirement that 
a <C d for the surface corrugations be small in height, 
used in the Sec. HVI when deriving the explicit form (|3.7p 
for the surface scattering potential "E,(k x , k' x ). 



VI. ATTENUATION LENGTH ANALYSIS 

In this section we discuss the attenuation length L n for 
three configurations of the waveguide. In the first part of 
the subsections we give explicit expressions for the 1 / L n 
for any form of the binary correlator W(x) . In the second 
part, we assume that the random surface profile £(x) has 
the Gaussian binary correlator, 

W(x) = cxp(-x 2 /2R 2 c ). (6.1) 

Then the RHP spectrum (|2.7p is given by 

W(k x ) = V2^R C exp(-k*R 2 c /2). (6.2) 

The RSGP spectrum T(k x ) defined by Eq. gl|), can be 
written as 



T{k x ) 



[{k x R c f - 4(k x R c ) 2 + 12] 



16i?3 
exp [-{k x R c ) 2 /A\ 



In what follows, we express the inverse attenuation 
length in the form of the dimcnsionless quantity A„/2L„. 
Since L n depends on as many as four dimensionless pa- 
rameters (ka) 2 , kR c , kd/ir, and n, a complete analysis 
appears to be quite complicated. For this reason, be- 
low we restrict ourselves to an analysis of the interplay 
between different scattering mechanisms, as a function 
of the dimensionless correlation length kR c for different 
values of (ka) 2 . We concentrate our attention on a mul- 
timode waveguide with a large number of propagating 
modes, N4 « kd/ir 3> 1. Our main interest is in two sit- 
uations, namely, in a small-scale roughness of the "white 
noise" kind for which kR c <C 1, and in large-scale random 
corrugations with kR c 3> 1. 



A. Waveguide with uncorrelated boundaries 



In view of Eq. (|4.20bp , the general expression (|5.2p for 
the inverse attenuation length shows that for a waveg- 
uide with uncorrelated boundaries this quantity consists 
of three terms, 



(1) 



L 



(2) 



'(2)1 



(6.4) 



These terms are originated from different mechanisms of 
the surface scattering. The first attenuation length iff 
is related to the RHP spectrum W(k x ) through the ex- 
pression for 2Q 1 (k x , k' x ). In accordance with Eq. (|4.17b|) 
it is given by 



1 

l-tn 



2a 2 (™/d) 2 g (,T,//,/l- 



k n d 



kn' d 



n' = \ 

x [W(k n + k n ,) + W(k n -k n ,)}. 



(6.5) 



Here the term 1/ Ln} n corresponding to n' — n, is related 
to the amplitude scattering while the terms , n 
result from the gradient scattering. These two mecha- 
nisms of surface scattering are due to the corresponding 
terms in the expression for Eu(k x , k' x ), see Eqs. (|4.11ap 
and (|4.12ap . The first one is related to the term depend- 
ing on the amplitude of the roughness profile a^i (x) , and 
the second one is due to the terms depending on the 
roughness gradient a^x). The expression (|6.5p exactly 
coincides with that previously obtained by various meth- 
ods (see, e.g., Ref. |l|). 

(2) 

The second attenuation length L„ ' related to the 
RSGP spectrum through 2Q2(k x ,k' x ), is associated solely 
with the SGS mechanism due to the operator Vi(x), see 
Eqs. (|4.13ap and (|4.14a|) . In accordance with Eq. (|4.18b[) . 
it is described by 



N d 



(6.3) 



(2) 



r (2) 
n' = l lj n,n' 



(6.6) 
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Here the terms in the sum are 



(im/d) 4 



r(2) 



(27m) 5 



[T(2k n ) + T(0)] 

(6.7) 



and 



4 16 (nn/df (nn'/d) 2 (n 2 + n' 2 ) 2 



'(2) 



(n 2 



(6.8) 



x [T(fc„ + fc n 0+T(fe n -fc„/)] 



The third term L„ 11 is also related to the RSGP spec- 
trum through Q-2\{k x , k' x ) due to the produc t £'^(x) ^(x), 
see Eq. (|4.15ap . In accordance with Eqs. (|4.19bj) . it is 
described by 



r(2)I 



E 



— r(2)I' 
n' = l ^njnf 



(6.9) 



where 



r (2)1 



a 4 (nn/d) 4 

: T kl 



2(nnf 



(6.10) 



and 



1 



L 



(2)1 



x [T(2fc„) + T(0)] . 

4 16 (vm/d) 2 (Tm'/rf) 2 (n 2 + n' 2 ) 2 
7T 4 fe„ fc n / (n 2 — n' 2 ) 4 



x [T(fc„ + fc„0 + T(fc„ - k„>)} cos 4 [7r(n - n')/2] 



.11) 



Analogously, the terms in Eq. (|6.4p can be classified as 
partial terms describing the intramode (1/Ln,li, 1/L^ n 
and 1/L^ n l ), and intermode {l/L ( ^ n ,_, n , l/£^ ¥n and 

(2) T 

V^nnVn) scattering. According to the discussion of 
Eq. (|6.5p , the amplitude and gradient terms result in the 
intramode and intermode scattering, respectively. As for 
the square-gradient terms, they lead to both types of 
scattering. 



1. Small-scale roughness 

To continue our analysis, let us start with the widely 
used case of a small-scale boundary perturbation, kR c <C 
1. In this case the surface roughness can be regarded 
as a delta-correlated random process with the correlator 
W(x — x') » W(Q)S(x — x') and constant power spec- 
trum W(k x ) as W(0) ~ R c . Taking into account the 
evident relationship fcA„ > 1, one can get the following 
inequalities specifying this case, 



kR c < 1 < kA r 



(6.12) 



It is necessary to underline that in the regime of small- 
scale roughness (|6.12|) the second of the weak-scattering 
conditions in Eq. ()5.3p is not so restrictive as the first one, 
and directly stems from it, R c -C A„ -C 2L n . 

In Eqs. f63|l . ([6J| . (|HTH|> . ([630]) and (fBTTTj) for the at- 
tenuation lengths the argument of the correlators W(k x ) 
and T{k x ) turns out to be much less than the scale of 
their decrease Rj 1 under the conditions (16.12|) . Thus, 
for any term in the summation over n' the argument can 
be taken as zero. 

Therefore, the first attenuation length is determined 
as follows, 



2L n 



, 2 n W(0) ^ (nn'/d) 2 
1 ' kd/ir k ^ k n ,d 



(ka) 



n kW(Q) 



kd/n 



(6.13a) 



(6.13b) 



The approximation is made here for a large number of 
the conducting modes, with the change of the summation 
over n' by integration. For a further estimate, one can 
take into account the formula 

/oo /*oo 
dxW{x) = 2R C / dpW(R c p), (6.14) 
-oo JO 

that directly follows from the definition (|2.7[) for the 
Fourier transform W(k x ) of the binary correlator W(x). 
The function W(R c p) is the dimensionless correlator of 
the dimensionless variable p, with the scale of decrease 
of the order of 1. As a result, the function W(R c p) does 
not depend on R c . Therefore, the integral over p enter- 
ing Eq. (|6.14p is a positive constant of the order of unity. 
For example, in the case of the Gaussian correlations (see 
Eq. ifO]) ). we have, W(0) = V2nR c , therefore, the inte- 
gral is \pK~j2. 

For the second length in Eq. (|6 .4[) one gets, 



A, 



Xka) 



T(0) 



2L^ 



k n d (kd/n) k 3 



16k n d 




(nn'/d) 2 (n 2 + n 12 ) 2 
k n ,d (n 2 - n' 2 ) 4 

n 2 (ka) 4 n 3 T(0) 
T k n d (kd/n) k 3 ' 



(6.15a) 



(6.15b) 



In Eq. (|6.15ap every term in the sum rapidly decreases 
with an increase of the absolute value of An = n — n' . 
This can be seen by making use of the following estimate, 



2 \2 



1 



(n 2 +n' ) 

(n 2 - n' 2 ) 4 ~ 4(An) 4 



for n»|An|. (6.16) 
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Therefore, the sum in Eq. (|6.15a|) can be correctly esti- 
mated by three terms with n' = n, n ± 1. For simplicity, 
in Eq. (|6.15b|) we assume Nd » n » 1, and replace the 
curly braces by factor 1/4. 

The explicit form for T(0) directly follows from the 
definition (14. 6| for the correlator T(k x ), 



no) 



dxW" 0) 



R.. 



dp 



d 2 W(R cP ) 
dp 2 



(6.17) 



If the roughness correlations are of the Gaussian form, 
then according to Eq. ([OD, we have T(0) = 3^/tt/4RI 
and the integral over p entering Eq. (|6.17p is equal to 
3^/4. 

The second term in Eq. (|6.4[) is bigger than the third 
one. Therefore, a good aproximation for Eq. (|6.4p is 



1 



L 



1 



(6.18) 



One can make sure of this if we consider in Eq. (|6.1ip 
the product of the fast decreasing factor (|6.16p and 
cos 4 [7r(n — n')/2] (the latter cancels the terms with An = 
1,3,..). Thus, we arrive at the following estimate, 



21 



(2)1 



An 
or (2)1 



(ka) 4 



no) 



36 k n d (kd/w) fc 3 



(6.19) 




FIG. 1: Plot of A n /2L ( n ) (solid line) and A n /2L%> + 



A n /2L<h )l (dashed line) versus kR c , for kd/n — 100.5, n — 50 
and (ka) 2 = 10" 2 . 



is due to the first attenuation length, Ln *C L)?' . The 
expression (j6.21[) shows that the crossing point is smaller 
for smaller values of ka, as well as for smaller mode in- 
dices n, or for larger values of the mode parameter kd/ir. 



-(2) 



By comparing Eq. (]6 . 19[) with Eq. (|6.15b[) one can con- 
firm Eq. (638]) ■ 

The expression (j6.18[) works well in the whole do- 
main of applicability of our theory, not only in the 
region of small-scale roughness. To support this, in 
Fig. Q] we compare the plot of A n /2L^ with the plot 

(21 (21T 

of A n /2L n + A n /2L n , they are practically the same. 
Therefore, in what follows for this case of uncorrelated 
boundaries we will neglect the third term in Eq. (|6.4p . 

Due to Eq. (|6 . . the relationship between the first 
and second attenuation lengths is expressed by 



L (1) n 2 



(6.20) 



According to a substantially different behavior of the 

quantities A n /2Ln^ and A n /2Ln' > with respect to kR c , 
it becomes clear that they have to intersect at the cross- 
ing point (kR c ) cr . If the crossing point falls onto the 
region of a small-scale roughness (kR c -C 1), its depen- 
dence on the model parameters is obtained by equating 
the expression (|6.20|) to unity, 



(kR c f cr ~ [ka)nh 



(6.21) 



To the left from this point (kR c ) cr the SGS length pre- 
vails, L { n ] < L ( n\ To its right the main contribution 



■scale roughness: weak correlations 



The intermediate situation arises when the correlation 
length R c becomes much larger than the wave length 
2-7r/fc, but still remains much less than the cycle length 
A„, 



(6.22) 



1 < kR c < kA, 



As before, the first of the weak-scattering conditions (|5.3[) 
is the most restrictive, R c -C A„ -C 2L n . 

Since the distance A„ is larger than the correlation 
length R c , the successive reflections of the waves from 
the rough surface are weakly correlated. Note that the 
distance between neighboring wave numbers k n and k n ±i 
is much smaller than the variation scale i?" 1 of the cor- 
relators W(k x ) and T(k x ), 



kn 



k n ±i \ « \dk n /dn\ = irAj < R~ 



(6.23) 



This implies that the correlators W(k x ) and T(k x ) are 
smooth functions of the summation index n' . Therefore, 
the sum in the expression (|6.5[) for the first attenuation 
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length Ln^ can be substituted by the integral, 



~ (ka^ 



(i) 



kdjix J 



2L,_ 

x\[W{k n + k n ,) + W(k n -k n ,)) 
k 



(kaf 



kdj-n J_ k 



dk x yk? 



W(k n - k x ) 
k 



(6.24a) 



(6.24b) 



Eq. (I6.24p shows that the first attenuation length is con- 
tributed by the scattering of a given n-th mode into other 
propagating modes. In order to obtain this asymptotic 
result, we have used only the condition of week correla- 
tions, R c <C A n . Therefore, Eq. (|6.24[) provides a reduc- 
tion to Eq. (|6 . 1 3[) for small-scale corrugations, kR c <C 1. 
In the case of a large-scale roughness, kR c ^> 1, the for- 
mula (|6.24p obtained for Ln allows one for further sim- 
plifications as was done in Ref. Q]. 

In contrast with L% , due to a rapidly decaying fac- 

________ 

tor (|6.16p the SGS length L n can be still described by 
keeping tree terms only, n' = n, n ± 1, in the sum in 
Eq. (|6.6p . Taking into account the estimate (|6 . 2311 . for 

the case Nd ^> n ^S> 1 one can obtain, 





_ 7T 2 {kerf 


n 3 


T(Q)+T(2A; n ) 






) '~ 2 k n d 


(kd/n) 


fc 3 




•{ 


'1 1 


2 8 
+ — 


1 


(6.25a) 


3 + (2t™) 2 _ 






_ 7T 2 (kaf n 3 


T(0)- 


1- T(2fc„) 


(6.25b) 
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To get the final expression (|6.25b|) . we replaced the curly 
braces in Eq. (|6.25ap by the factor 1/4. Naturally, at 
small-scale corrugations, kR c <C 1, the obtained result 
(|6.25p reduces to Eq. (|6.15[) . For a large-scale roughness, 
kR c 3> 1, one should use Eq. (|6.25p because of an arbi- 
trary value of the parameter k n R c . We do not consider 
here this case in detail due to its intermediate character. 



3. Large-scale roughness: strong correlations 

Another limit case refers to the correlation length R c 
to be very large both in comparison with the wave length 
27r/fc and with the cycle length A„, 



1 < kK n < kR c . 



(6.26) 



In this case the number of wave reflections on the scale 
of the correlation length R c is large. Therefore, the suc- 
cessive reflections are strongly correlated to each other. 

Under the relations (|6.26p the second of the weak- 
scattering conditions (|5.3p is the most restrictive, there- 
fore, the condition of the applicability reads as 



A„ < R c < 2L r . 



(6.27) 



The latter requirement determines the upper limit for the 
value of the correlation length R c . 

Due to Eq. (|6.26|1 . the distance between neighboring 
wave numbers k n and k n ±i turns out to be much larger 
than the variation scale R~ x of the correlators W(k x ) 
and T(kx), 

\k n - k n±1 \ « \8k n /dn\ = ttA- 1 » R-\ (6.28) 

This indicates that the probability of the intermode 
(p! 7^ n) transitions is exponentially small and the at- 
tenuation lengths are mainly formed by the incoherent 
intramode (n' — n) scattering. Formally, for strong 
correlations (|6.26p the correlators W(k x ) and T{k x ) are 
sharpest functions of the summation index n . In the 
sums of Eqs. (|6.5p and (16. 6p for the attenuation lengths 
the main contribution is due to the diagonal terms with 
n' = n, for which W(2k n ) and T(2k n ) can be neglected 
in comparison with W(0) and T(0). 

Thus, the first attenuation length reads 



A r 



{kaf 



kW(0). 



2L W k n d (fcd/vr) 3 
Correspondingly, for the SGS length one gets, 
A„ it 2 (kaf n 3 T(0) " 



2L 



(2) 



2 k n d (kd/%) k 3 



1 



(2nnf 



(6.29) 



(6.30) 



The ratio of the first attenuation length to the second 
one can be presented as 



r(l) 

7w 



A, 



2L 



(i) 



(k n K r 



< 1. 



.31) 



According to the first inequality in Eq. (|5.3p . applied to 
the condition (|6 . 26(1 of the strong-correlations and to the 
evident relationship k n h. n > 1, one can see that the am- 
plitude scattering length always prevails over the SGS 
length within the interval of strong correlations. For this 
reason in the condition of applicability (|6.27p one should 

substitute L n by Ln ■ This allows one to arrive at the 
inequalities in the explicit form, 

kA n < kR c < (kknfka^ikd/im). (6.32) 

Before we pass to the next configuration of the waveg- 
uide, it is interesting to compare the case of two uncor- 
rected boundaries with the waveguide with one rough 
boundary, studied in Ref. [17[. One can found that un- 



der the approximation (|6.18p . the expression for 1/L r 
turns out to be as twice as bigger than that one for the 
one-rough-boundary configuration. This fact is quite nat- 
ural since in the case of uncorrelated boundaries the two 
rough surfaces are of the same statistical nature (see the 
factor 2 in the first two terms of Eq. (|4.20bp ). 

The remarkable similarity between the behavior of the 
inverse attenuation length for both configurations can be 
observed in Fig. [21 Note that the normalization factor A n 
for one boundary is twice as much of the quantity for two 
boundaries, and both curves are practically superposed. 
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FIG. 2: Plot of A„/2L„ versus kR c for kd/n = 100.5, n = 50 
and different values of (ka) 2 . Different sets of curves are 
shown for three waveguide configurations, as well as for the 
waveguide with one rough boundary analyzed in Ref. fl7l ] 
(dashed line). Curves without symbols depict Eq. (16. 4[) for 
the waveguide with uncorrelated boundaries; curves labelled 
by triangles and diamonds stand for Eqs. (|6.37[) and (|6.44[) for 
waveguides with antisymmetric and symmetric boundaries, 
respectively 



modes for which the difference of the indexes n — n' is 
an even number. There are transitions only between odd 
and even channels. This fact can be explained by an ex- 
istence of an additional integral of motion, namely, the 
parity. Specifically, if one solves the stationary problem 
for the effective Hamiltonian, there are two sets of eigen- 
states of different parity. Correspondingly, there are two 
independent sets of eigenvalues, see Refs. (l6l. [l8|. 

The intramode scattering, related to the RSGP spec- 
trum through the kernel Qi{k x ,q x ) in the Eq. (|4.27bj) . is 
related to the SGS mechanism via the function V(x) in 
Eq. (|4.24ap . Therefore, we have 



1 



(nn/d) 4 



[T(2k n ) + T(0)] . 



(6.35) 



Note that in a single-mode waveguide with = 1 the 
sum over n' in Eq. (|6.33p contains only one term with 

n' = n = 1. In this case the backscattering length 
enterin g in Eq. (|6.35|) is in accordance with that obtained 
in Ref. [II. 

Bearing in mind very different dependencies of the 
W(k x ) and T(k x ) power spectra on the roughness cor- 
relation length i? c , one can write the inverse attenuation 
length as the sum of two terms. One is related to W(k x ), 



B. 



Waveguide with antisymmetric rough 
boundaries 



According to Eqs. (|4.28b|) and (pT2"]) . the inverse atten- 
uation length can be presented as a sum of the partial 
inverse lengths \jL n . n i that are responsible for the scat- 
tering from rtth into n'th modes, 



1 



E 



(6.33) 



However, it is noteworthy that in this case the partial 
terms describing the intramode \jL n ^ n and intermode 
scattering are rather different since they are 
associated with different scattering mechanisms. Indeed, 
the intermode scattering results from the RHP spectrum 
only, through the expression in Eq. (|4.26b[) for the ker- 
nel Qi(k x , q x ). These partial lengths are formed by the 
gradient scattering mechanism, due to the term contain- 
ing in the scattering potential the derivative £,'(x), see 
Eq. (14.23a|) . They read as 



1 



4cr 



(irn/d) 2 {irn'/d) 2 



fan.' d 



[ir[n-n')/2] 



x[W{k n + k n ,) + W(k n -k n ,)}. (6.34) 

From the above expression one can draw an interest- 
ing conclusion. As one can see, in the case of symmet- 
ric boundaries the transitions are forbidden between the 



1 



Nd 

E 

n' = 1, 
n' =/= n 



1 



(6.36) 



and the other is just the inverse of the intramode at- 
tenuation length determined by Eq. (|6.35|) . Therefore, 
1/L n = 1/L nj „, and one can write, 



1 



L 



(1) 



r(2)- 



(6.37) 



In Fig. [5] we depict the length defined by Eq. (|6.37[) as 
a function of kR c for different values of ka; the curve is 
labelled with triangles. 

As above, it is convenient to consider separately three 
regions of the correlation length. 



1. Small-scale roughness 

We start with the standard case of a small-scale bound- 
ary perturbation specified by the conditions (|6.12p . In 
this case for any term in the summation over n' the ar- 
gument of the correlators W(k x ) and T(k x ) is much less 
than the scale of their decrease Rj 1 , therefore, it can be 
taken as zero. Thus, for the attenuation length describ- 
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ing the intermode scattering one gets, 
A„ o n W(0) 



2L 



(i) 



A(ka) 



kd/n k 



g (nn'/df 

'-1 



n 



(kaf 



kW(0) 



kd/n 



[n(n-n')/2] (6.38a) 



(6.38b) 



The factor sin 4 [7r(n — n')/2] gives alternating values of 
one and zero for n — n' = 1,3, .. and n — n' = 2, 4, .., 
respectively, therefore, for a large number of conducting 
modes Nd ^ 1 it can be replaced by 1/2. The estimate 
()6.38bp is obtained with the change of the summation by 
integration. For a further analysis one has to take into 
account the formula in Eq. (|6.14|) for VF(0). 
The second term in Eq. (|6.37|) reads as 



A n n 2 {kaf 



T(0) 



2L n 



2 k n d (kd/n) k 3 



(6.39) 



where the explicit form for T(0) can be taken as in 
Eq. f6T7| . 

It is easy to realize that the relationship between 
lengths defined by Eqs. (|6.38b[) and (jQ9)) is of the 
same order as in the case of uncorrelated boundaries, 
see Eq. (|6.20[) . In this way one obtains the expression 
(|6.2ip for the crossing point (kR c ) cr at which the curves 
of A n /2Ln' > and A n /2L^ intercept each other. 



2. Large-scale roughness: weak correlations 



This intermediate situation is specified by Eq. (|6.22[> . 
In Eq. (|6.36p we replace sin [n (n — n')/2] by 1/2, and 
change the sum by the integral since the correlator W(k x ) 
is a smooth function of n' (see the discussion before 
Eq. (I6» . 



A,, 



(ha) 



2L 



(i) 



kd/n 



W(k n - k x ) 



(6.40) 

Note that Eq. (jUH)]) coincides with Eq. (|6.24bj) for the 
case of two uncorrelated boundaries. 

(2) 

For the length L n in the large-scale-roughness region, 
kR c 3> 1, one should use the general expression (|6.35|) 
because of arbitrary value of the parameter k n R c . Its 
normalized version reads 



A, 



n 2 (kaf n 3 T(2k n ) + T(0) 



2L ntn 4 k n d kd/n 



k 3 



(6.41) 



3. Large-scale roughness: strong correlations 



In this limit case specified by Eq. (|6.26[) . the distance 
between wave numbers k n and k n ±\ is much larger than 



the variation scale R~ x of the correlator W(k x ), see 
Eq. (|6.28p . This means that the probability of the in- 
termode (n 1 7^ n) transitions is exponentially small and 
the total attenuation length is mainly formed by the in- 
coherent intramode (n' = n) scattering, 



1 



1 

-12)"' 



(6.42) 



For strong correlations the correlator T(k x ) is a sharp 
function of its argument k x . Thus, the term with T(2k n ) 
can be neglected in comparison with T(0), and one gets, 



A„ 
2L n 



(ka) 



T(0) 



4 k n d (kd/n) k 3 



(6.43) 



In this region of kR c the SGS mechanism is recovered, 
determining the total attenuation length. The second 
inequality in Eq. (|6.27p is automatically satisfied for any 
value of kR c 



C. Waveguide with Symmetric Rough Boundaries 

As above, for this configuration it is convenient to write 
down the inverse attenuation lengths as a sum of two 
terms, 



1 



i 



l 



(6.44) 



Here each term is associated with the corresponding term 
in Eq. (|4.36b|) . resulting from different mechanisms of the 
surface scattering. 

The first term in Eq. (|6.44[) is related to the RHP spec- 
trum W(k x ) through the expression for Qi(k x ,k' x ). Ac- 
cording to Eq. (|4.34bp , we have 



-(i) 



= 4cr 



2 (nn/d) 2 
kn.d 



E 



(nn'/d) 2 
k 71 /d 



cos 



[n(n - ri) /2] 



x [W(k n + Av) + W(k n - h*,)] 



(6.45) 



The term 1/Ln}n corresponding to n' = n, is related to 



the amplitude scattering while the terms 1 /L^' n ,, n result 
from the gradient scattering. These two mechanisms of 
surface scattering are due to the corresponding terms in 
the expression for Ei(k x , k' x ), see Eq. (I4.32p . Namely, the 
former follows from the term depending on the amplitude 
of the roughness profile er£j (x) , and the latter from the 
terms depending on the roughness gradient cr^(x). 

The second term in Eq. (jolgp related to the RSGP 
spectrum through Qi(k Xl k' x ), is associated solely with 
the SGS mechanism due to the square-gradient function 
V(x), see Eq. (14331) . In accordance with Eq. (|4.35b|) . it 
is described by 



(2) 



r (2) 
n' = l lj n,n' 



(6.46) 
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Here 



-(2) 



a 4 {im/df 
T kl 



(7m) 5 



[T(2k n ) + T(0)] 

(6.47) 



and 



128a 4 (nri/d) 2 (rrn'/d) 2 (n 2 + n l2 f 



-(2) 



k n > {n 2 - n' 2 Y 



x cos 4 [7r(n — n')/2] 

x [T(fe„ + k n ,) + T(k n - k n ,)] 



(6.48) 



Note that in a single-mode waveguide with Nd = 1 
(therefore with n' = n = 1 in Eqs. (|6.45[) and (|6.46|) ) 
the backscattering length takes the expression de- 
rived in Ref. (l5j . 

The diagonal part of the inverse attenuation length, 



l/£n,n = l/i^.n + l/in%, characterizes the intramode 
scattering. The off-diagonal partial attenuation length 



-(2) 



scattering 



l/L 



(2) 



describes the intermode 



To compare with the above case of the antisymmet- 
ric boundaries, here there are no transitions between the 
modes with odd values for the difference n — n' . There- 
fore, only the transitions between odd or even modes are 
allowed only. This fact is related to the underlying sym- 
metry, and can be associated with specific integrals of 
motion of the total Hamiltonian. 



1. Small-scale roughness 



In the case of a small-scale roughness the first attenu- 
ation length reads 



2L« ~ [ ' kd/n k 2^ k n ,d 
n fcW(O) 
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For the second length we have, 
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(6.50a) 



(6.50b) 



As well as in the configuration with uncorrelated bound- 
aries, due to the estimate (|6.16[) every term in the sum 
of Eq. (|6.50a[) rapidly decreases with an increase of the 
absolute value of Art = n — n' . In addition, the factor 
cos 4 7r(n — n')/2 is zero for n — n' = 1,3, ... Therefore, the 
sum in Eq. (|6.50a|) can be evaluated by three terms with 
n' = n, n ± 2, and one can replace the quantity in the 
curly braces by 1/4. Here we also assume, Nd ^> n ^> 1. 

It is evident that in the small-roughness case the cross- 
ing point is determined by Eq. (|6.21|) . 



2. Large-scale roughness: weak correlations 



In the intermediate situation (|6.22|) the first attenu- 



ation length Ln coincides with that obtained for the 

W(k n - k x ) 



previous configurations, 
A„ (fee) 2 n 
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(2) 

As for the second length L n , for Nd n 3> 1 it can be 
written, 
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Tx 2 (kof n 3 T(0) + T{2k n ) 
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3. Large-scale roughness: strong correlations 

In this region of kR c the first attenuation length reads, 
A„ Aka) 2 n 3 



2L W k n d (kd/TT) 3 

For the SGS length one gets, 
A n vr 2 ^) 4 n 3 T(0) 
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The ratio of the first attenuation length to the sec- 
ond one can be obtained due to Eq. (|6.31[) . One can 
see that the amplitude scattering length always prevails 
over the SGS length within the interval of strong corre- 
lations, L [ n ] /L ( n ) < 1. For this reason, the condition of 
applicability coincides with Eq. ()6.32|) . 

In Fig.[2]by diamonds we show the behavior of the total 
dimensionless inverse scattering length as a function of 
the correlation parameter kR c . 



VII. NUMERICAL ANALYSIS OF THE 
SURFACE SCATTERING POTENTIAL U 

In this Section we present numerical data for periodic 
waveguides (billiards) with the potential U(x,z). Our 
main attention is payed to the structure of eigenstates for 
the discussed types of boundaries. Although this infor- 
mation is not directly related to the scattering properties 
of open waveguides, it is instructive to reveal a role of the 
SGS-terms in the potential, by comparing the structure 
of eigenstates for different boundaries. The numerical 
method is described in details in Ref.[l3|; below we dis- 
cuss the most interesting results and relate them to the 
previous studies (see, Ref.[i~?| and references therein). 

Our particular interest is in the "repulsion" effect dis- 
covered in Ref. [l8[ for quasi-lD billiards with rough sur- 
faces. This phenomena can be observed for a subset of 
eigenfunctions with the smallest values of the transver- 
sal quantum number of the unperturbed billiard with 
perfectly flat boundaries. For these eigenstates it was 
found that the maximum of the eigenfunction intensity is 
strongly shifted towards from the rough boundary. This 
"repulsion" was shown to increase with a decrease of the 
roughness correlation length R c . Based on our numerical 
data, one can make a conclusion that this effect of "re- 
pulsion" is directly related to the SGS-terms discussed 
above. 

To start with, we show in Fig. [3] a typical eigenstate of 
this type for the billiard with one rough surface (upper 
panel). To compare with, the low panel presents the 
corresponding eigenstate for the same potential, however, 
without SGS-terms. As one can see, by removing these 
terms, the repulsion disappears. 

It is also instructive to compare the structure of such 
eigenstates with those computed for the case of both sym- 
metric and antisymmetric rough boundaries, see Fig. [4] 
In both cases the potentials have the SGS-terms. How- 
ever, the repulsion phenomena is absent in the antisym- 
metric billiard. A close inspection reveals that the dif- 
ference between these two cases is entirely due to the 
fact that in the antisymmetric billiard there is no square- 
gradient scattering between different channels, in con- 
trast to the symmetric billiard. Thus one can make a con- 
clusion that the main contribution to the repulsion is due 
to the intermode SGS-terms. This conclusion supports 
an observation found in Ref. [18| according to which the 
effect of the repulsion can be explained due to a strong lo- 
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FIG. 3: (Color online) Upper panel: intensity of a typical 
eigenstate manifesting a "repulsion" from the rough surface. 
Low panel: the corresponding eigenstate for the potential 
without the SGS-terms. 



calization in the channel space. This localization occurs 
due to a relatively strong interaction between different 
conducting channels (waveguide modes). 



VIII. CONCLUSION 

We have developed a perturbative approach to the sur- 
face scattering in multimode quasi- ID waveguides with 
rough surfaces. The approach is based on the con- 
struction of an effective Hamiltonian that formally corre- 
sponds to the "bulk" scattering described by the poten- 
tial of a very complicated structure. The detailed analy- 
sis of this potential has led to the conclusion that there 
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FIG. 4: (Color online) Intensity of an eigenfunction reveal- 
ing strong repulsion from two symmetric rough boundaries 
(low panel). The same for the antisymmetric billiard for 
which there is no interaction between different channels (up- 
per panel). 



are specific square-gradient scattering (SGS) terms that 
were not taken into account in the previous studies. We 
show that in spite of their seeming smallness in com- 
parison with the standard terms, they can give a major 
contribution to the value of the scattering length in the 
case of a sufficiently strong roughness characterized by 
a corresponding value of the derivative of boundary pro- 
files. 

In order to demonstrate a strong influence of this SGS- 
mechanism, we have applied the approach to the waveg- 
uides with the boundaries of three different types. Specif- 
ically, we considered the waveguides with two uncorre- 
cted boundaries, and compared the results with those 
obtained for other two cases of symmetric and antisym- 



metric rough boundaries. In this way one can see the 
influence of various terms very clearly. 

According to our analysis, all terms in the expression 
for the effective potential can be regrouped in three parts. 
The first one is due to the roughness mean amplitude a. 
Another group is characterized by the derivative of pro- 
files, reflecting the fact that the scattering depends on 
two characteristic parameters, namely, on the amplitude 
a (or, the variance a 2 ), and on the correlation length R c 
inversely proportional to the mean derivative of the pro- 
files. In the expression for the scattering potential the 
above two types of terms have linear dependence on the 
root-mean-square roughness height a. There is one more 
group of terms that can not be reduced to the two previ- 
ous ones, and that gives rise to the square dependence on 
the profile derivative a/R c . For this reason, in the previ- 
ous studies it was assumed that these SGS-terms can be 
safely neglected. 

However, we show that apart from the cr-dependence 
due to the amplitude and gradient terms, as well as the 
cx 2 -dependence due to the SGS-terms, there is also depen- 
dence on the correlation parameter R c . As was found, the 
i? c -dependence for the SGS-terms turns out to be very 
different from those corresponding to the amplitude and 
gradient terms. Our analysis reveals that there is a region 
of parameters where the SGS-contribution to the inverse 
scattering length is compared (or even larger) with the 
standard contributions due to the amplitude and gradi- 
ent scattering. 

In accordance with the general theory of surface scat- 
tering [l[ the wave propagation through any n-th con- 
ducting channel is determined by its total (or outgoing) 
attenuation length L n . The inverse attenuation length 
l/L n can be presented as the sum of partial inverse 
lengths, see Eqs. (|rT3j) . (|OTj) and (HU). The first term 

is associated with both the amplitude and gradi- 
ent scattering, it depends on the roughness-height power 
spectrum W(k x ). To the contrary, the partial length 

(2) 

L n being connected with the SGS-terms, depends on 
the power spectrum T(k x ) that is principally different 
from W{k x ). For the convenience of comparison between 
three different types of boundaries, in Eq. (|6.4p we have 

(2) T 

separately written the term 1/L n that appears in the 
specific case of two uncorrelated boundaries. This term 
also depends on T{k x ). 

Such a representation allows us to analyze the total 
dependence of the inverse attenuation length l/L n on 
two control parameters, the roughness hight a and the 
correlation length R c . Our analysis shows that one can 
speak about three characteristic regions of R Cl namely, 
the small-scale (16. 12)) . large-scale weak-correlation (|6.22|) 
and large-scale strong-correlation (|6.26[) roughness. 

As we noticed above, the W(k x )- and T(fca;)-power 
spectra have very different dependencies on the rough- 
ness correlation length R c . Specifically, the inverse value 
of the L2-length decreases with an increase of the corre- 
lation parameter kR c ; this occurs in the whole interval 
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of the roughness values. As for the inverse value of the 
Li-length, it shows a different behavior in each of the 
three regions of roughness and strongly depends on the 
waveguide configuration. 

In the region of a small-scale roughness (|6.12[) the 
inverse length 1/Li increases with an increase of kR c . 
Therefore, the curves corresponding to the terms l/Ln 

and (l/Ln + l/Ln for the case of uncorre- 

cted boundaries) intersect each other, see Fig. 2. If the 
crossing point (kR c ) cr falls into this region of small-scale 
roughness, it obeys the law given by Eq. (|6.21| . To the 
left from the crossing point (kR c ) cr , the SGS-length L 2 
prevails over the standard length, Lffl <C L„ ■ Con- 
trary, to the right from the point (kR c ) cr , the Li-length 
gives the main contribution to the scattering process, 

Lit' < L 

The detailed analysis of three types of boundaries 
shows that for the small-scale roughness region one can 
specify that the crossover (with an increase of kR c ) oc- 
curs due to the following interplays: 

• uncorrelated boundaries: the two SGS-terms that 
depend on the T(k x )-power spectrum, versus the 
amplitude and gradient terms depending on the 
W (k x )-powev spectrum; 

• antisymmetric boundaries: the SGS-term versus 
the gradient term; 

• symmetric boundaries: the SGS-term versus the 
amplitude and gradient terms. 

Thus, we show that in the region of small-scale rough- 
ness and for any fixed value of the roughness height a, 
one can indicate the region of small values of the correla- 

(2) 

tion length R c where the SGS-length L n predominates 
over the the length lJ$ . This predominance occurs in 
spite of the fact that 1/Zn is proportional to a 1 while 

1/Ln is proportional to a 4 . 

It is interesting to note that in the region to the 
left from the point (kR c ) cr we have observed that the 
length L^ B for symmetric boundaries is twice as much 
as the length L^ B for the uncorrelated boundaries. On 
the other hand, the latter length is twice as much as 
the length L^ B for the antisymmetric boundaries, i.e., 



There is another border that is characterized by the 
transition to the large-scale roughness with strong cor- 
relations, where a very different behavior of the total 
inverse length 1/L n was found for the antisymmetric 
boundaries, in comparison with the uncorrelated and 
symmetric ones. For the two last cases, the major mech- 
anism of the scattering is due to the amplitude terms in 
the potential, and a similar behavior with respect to kR c 
(l/L n oc kR c ) occurs for both of them. However, the 
attenuation length for the uncorrelated boundaries turns 
out to be twice as much as the length for the symmetric 
boundaries, i.e., L^ B rj 2L^ b (note that this relation 
is inverted in comparison with the region of small-scale 
roughness). 

As for the waveguide with the antisymmetric bound- 
aries for which the amplitude terms are absent in the po- 
tential, the SGS-mechanism prevails in this region, lead- 
ing to the dependence 1/L n oc (kR c )~ 3 . Note, that here 
there is the interplay between the SGS-term and gradient 
one. The transition point was found to be kR c ~ kA n , 
with the crossover specified by the following interplays: 

• uncorrelated and symmetric boundaries: the am- 
plitude and gradient terms versus the amplitude 
terms, 

• antisymmetric boundaries: the gradient and SGS 
terms versus SGS-term. 

Finally, there is an intermediate region of a large-scale 
roughness and weak correlations where the first attenu- 
ation length mainly contributes, L« <C Ln . Here one 
can observe a strong similarity of the inverse attenua- 
tion length behavior for all configurations considered. It 
should be noted that this region is quite difficult for ob- 
taining simple analytical estimates due to a strong depen- 
dence of the number of terms in the corresponding sums 
on the parameter kR c . Numerical data manifest that the 
attenuation length roughly remains the same inside this 
intermediate region. 
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